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q Abstract 

We develop a new Voronoi protocol, which is a space tessellation method, to generate a fully dense (containing no voids) model 
of nanocrystalline copper with precise grain size control; we also perform uniaxial tensile tests using molecular dynamical (MD) 
simulations to measure the elastic moduli of the grain boundary and the grain interior components at 300 K. We find that the grain 
^< boundary deforms more locally compared with the grain core region under thermal vibrations and is elastically less stiff than the 
core component at finite temperature. The elastic modulus of the grain boundary is lower than 30% of that of the grain interior. Our 
' results will aid in the development of more accurate continuum models of nanocrystalline metals. 
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1. Introduction 

Solid metallic systems usually possess polycrystalline struc- 
tures composed of crystalline grains of different sizes and ori- 
entations. Polycrystalline metals with a grain size smaller than 
100 nanometers are called nanocrystalline metals. Especially, 
because of their small grain size, the grain boundary component 
occupies a substantial part in these materials. We can think of 
nanocrystalline metals as composites of crystalline grain cores 
! and grain boundaries. A grain boundary is the interface of finite 
average thickness and nonzero area between two neighbouring 
crystalline cores. Both the grain core and the grain boundary 
components play important roles in determining the bulk me- 
chanical properties of nanocrystalline metals [jl], |2|, |3|] . In this 
study, we focus on the elastic modulus of the grain boundary 
component of nanocrystalline metals. Understanding the elas- 
tic behaviour of the grain boundary region would be beneficial 
in evaluating the stress field around crack tips and dislocations, 
and can help elucidate the effect of porosity on the mechanical 
properties of nanocrystalline metals. 

In the past twenty years, it has been widely accepted that the 
grain boundary is about 70 to 75% as stiff as the grain core 
component J3,|5t]. However, this belief is based on numerical 
studies of specific grain boundaries, such as the relatively stable 
25 twist boundary, where atoms interact with one another via 
simplified Lennard- Jones potential at zero temperature , or 
on experimental studies assuming that the grain boundary com- 
ponent behaves like amorphous alloys |4j, |2t ISl I2D ■ A thorough 
investigation from a general atomic scale structure at finite tem- 
perature is still lacking. 
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In general, mechanical properties of nanocrystalline metals 
are affected by both grain size d lflo[[lTll and its dispersity cr s {d) 
1 1211 . Unfortunately, the microscopic structure of nanocrys- 
talline metals has not been fully characterized experimentally, 
and different fabrication processes alter the microstructure ex- 
tensively. Typically, Voronoi construction with randomly cho- 
sen Voronoi seeds is a model generating nonuniform grains and 
resembling nanocrystalline metals made by inert gas condensa- 
tion 11311 . Here, we use Voronoi seeds from random close pack- 
ing (RCP) (by that, we focus on its randomness and packing 
density of about 0.64, not its strict mathematical definition 1 14]) 
of identical spheres to construct a new model of fully dense 
nanocrystalline copper with uniform grains and a well-defined 
d, the diameter of a spherical grain approximating the polyhe- 
dral grain created by the algorithm, so that we can investigate 
the grain size dependence of the elastic modulus E separately, 
while ignoring the effect of cr s {d), a setup still experimentally 
unattainable, and has not been achieved by previous simula- 
tion studies. Using randomly packed uniform grains makes sure 
our model is not weakened on the granular scale due to the in- 
troduction of equal-sized grains and their possible ordered ar- 
rangement, and enables us to focus on the sub-granular struc- 
ture such as the grain boundary component. Knowing this is 
an ideal model, we carefully establish its validity by inspecting 
its microstructure such as grain boundary thickness, the overall 
elastic and plastic behaviors. The results of all these tests agree 
very well with what have been reported in simulations and ex- 
periments, and therefore we believe our model is adequate to 
analyze the elastic modulus of the grain boundary component, 
a quantity beyond the approach of simplified bicrystal model or 
full-sized experimental measurement. 

To study the elastic behaviour, we performed 3D simulations 
of uniaxial tensile tests at 300 K, about 22% of the melting tem- 
perature of copper. It is essential that this temperature is low, 
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so that our models do not alter their structures dramatically dur- 
ing tensile tests; on the other hand, it is high enough for us to 
observe the thermal effect on the elastic moduli, and there ex- 
ist many other experimental results at similar temperatures for 
comparison. We focused on systems with d smaller than 25 nm. 
Larger than this value, the grain core component occupies the 
major part of the system and the elastic contribution of the grain 
boundary component is negligible [l5t] . 






2. Nanocrystalline Model of Randomly Packed Uniform 
Grains and Molecular Dynamical Simulations 

The uniaxial tensile test has been considered the most di- 
rect way of determining the mechanical properties of a ma- 
terial 111 611 . We implemented a new Voronoi-like algorithm to 
generate the polycrystalline initial configurations for molecular 
dynamics (MD) tensile tests. Unlike the conventional Voronoi 
algorithm, where positions of Voronoi seeds are randomly cho- 
sen, leaving d ill-defined due to its high dispersity, particularly 
when the number of grains is small, the positions of seeds in 
our algorithm are mapped from center positions of 3D random 
close packings (RCPs) of monodisperse hard spheres, whose 
volume fraction is close to the random close packing density, 
* 0.64 m. Because the average distance between any pair of 
Voronoi seeds is nearly a constant in this algorithm, it allows us 
to generate almost evenly-sized Voronoi grains and control the 
dispersity of d accurately. 

2.1. Random Close Packing Finder of Identical Spheres 

First we generate random close packings of identical spheres 
under periodic boundary conditions. The algorithm begins with 
a non-overlapped random collection of spheres in a unit cell 
whose packing density is at least two orders of magnitude lower 
than the 3D random close packing density, ss 0.64. The inter- 
action between spheres can be described by the soft, repulsive 
potential 



V,(ry) 
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where e is the characteristic energy scale, r,-y is the distance 
between the z'th sphere of diameter Rj and the jth sphere of di- 
ameter Rj, and R,j = Rj + Rj. The average potential energy of 
the system is given by the expression 
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where n is the total number of spherical grains in the system. 
After an initial configuration is generated, we grow (shrink) 
each sphere by AR so that the volume fraction <p of the sys- 
tem increases (decreases) by 10~ 4 , if \ t0 , is smaller (greater) 
than a minimal average potential energy V m ,„ = 10~ 16 given by 
the machine precision. Each change of <p is followed by en- 
ergy minimization of V„„ done by the conjugate gradient (CG) 
method. During the CG relaxation process, spheres are allowed 
to move freely. AR/2 instead of AR is applied when switching 



Figure 1: Algorithm for random close packings of identical spherical particles 
shown schematically in 2D. Intermediate states without overlap (al) will grow, 
and overlapped intermediate states (a2) will shrink until a state (b) is reached. 
During the packing process, particles are allowed to move freely. 



from compression to expansion of the system and vice versa. 
The procedure of energy minimization following volume per- 
turbation is repeated until V„„„ < \, ot < 2V„„„ to create a pack- 
ing, where there are only tiny overlaps between spheres and 
each sphere keeps force balance with all its contact neighbors. 
The algorithm is schematically shown in figure [Till ST. 1 1911. 

The algorithm runs at zero temperature and under zero grav- 
ity, so the system is very unlikely to be trapped by global min- 
ima such as the fee or hep configuration on the potential energy 
landscape. In some rare cases, the packing algorithm ends up 
with unstable or ordered configurations, which can be removed 
manually. The number of grains n and volume fraction <p of 
the specific sphere packings used to obtain the Voronoi seeds 
in this study are (n,0) = (8,0.64), (27,0.63), and (512,0.64), 
respectively. All these sphere packings are randomly packed, 
and their <p's are basically the same as the RCP density. 



2.2. Algorithm for Voronoi Tessellation 

The centers of randomly close packed spheres are used as 
Voronoi seeds to create Voronoi grains by the standard Voronoi 
tessellation: the ;th Voronoi grain contains points that are closer 
to the /th seed than to any other seeds. The Voronoi boundary 
between the ;th and y'th spheres is determined by the relation 



(3) 



where r v ; is the distance between a point v in the Voronoi 
boundary to the ;th sphere of radius Rj. We call Voronoi seeds 
obtained this way 'uniform Voronoi seeds' and the created 
Voronoi grains 'uniform Voronoi grains'. In figure |2] it has 
been shown that the grain size dispersity <r s (d) ~ cr s (V vorono i) 
of uniform Voronoi grains is very small (< 5%), compared with 
that of randomly chosen nonuniform Voronoi grains (> 36%), 
and therefore influence of grain size dispersity is quantitatively 
negligible. Figure [3] shows the Voronoi tessellations of N - 
8,27,512 used specifically in this study. We fill each Voronoi 
grain with a face-centered cubic (fee) randomly-oriented crys- 
tallite. Using uniform Voronoi seeds guarantees that not a sin- 
gle Voronoi grain is considerably larger or smaller than the aver- 
age grain size, and each grain is surrounded by grain boundaries 
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Fi gure 2: Normalized standard deviation o~ s of Voronoi volume Vvoronoi 

for 

a single set of uniform Voronoi grains (black) and randomly chosen nonuni- 
form grains (blue) in a unit cell. For all values of n (8, 27, 512) used in this 
study, <r s (y vorono i)l V vorono i is smaller than 0.05 for uniform Voronoi grains, but 
greater than 0.36 for nonuniform grains. Volumes of grains are calculated using 
Voro++ (U]. 




Figure 4: A new Voronoi-like algorithm for generating nanocrystalline config- 
urations: (a) First, create a packing of monodisperse spheres. Its schematic 2D 
section is shown here. Positions of the spheres (orange dots) are used as seeds 
in the standard Voronoi procedure to create uniform Voronoi grains, (b) Second, 
fill the uniform Voronoi grains with randomly-oriented fee crystals (represented 
by black arrows), followed by MD relaxation at finite temperature. Each grain 
will then be surrounded by uniform grain boundaries (blue lines) with finite 
thickness. The size of the supercell has to be changed to get the desired grain 
size. We approximate these polyhedral Voronoi grains by spheres with an aver- 
age diameter d and an average grain boundary thickness cIqb- The overall grain 
boundary thickness of the system is 2dos- 




5 x 10 8 s _1 . Young's modulus is strain-rate independent though. 
We calculated the overall Young's modulus E using data in the 
strain (e) interval e < 0.3%, which is clearly within the linear 
stress(cr)-strain(e) region 111 311 . 




Figure 3: Voronoi tessellations of randomly close packed identical spheres of 
(a) n = 8, (b) n = 27, and (c) n = 512 used specifically in this study. Edges of 
Voronoi boundaries are shown by blue lines. Voronoi tessellations are visual- 
ized using Voro++ l2(ill . 



with similar geometrical features. In other words, the configu- 
ration can be treated as a homogeneous medium down to the 
characteristic length scale d in this study. 

Moreover, since the structure is randomly packed, it dis- 
misses the concern that there may exist any artificial slip planes 
in the system caused by orderly arranged grains with a uni- 
form grain size distribution, which may weaken the stiffness 
and strength of the system. After an initial configuration is con- 
structed, under zero external loadings, we relax it at the tem- 
perature that is also used for the subsequent tensile tests, until 
the MD extended hamiltonian reaches a stable value. This re- 
laxation process usually takes less than 1CT 10 seconds. The new 
Voronoi-like algorithm is schematically shown in figure |4] 

2.3. Molecular Dynamical Simulation of Tensile Tests 

Using the Parallel Molecular Dynamics Stencil (PMDS) 
code 1,21], we performed three dimensional MD uniaxial ten- 
sile tests, with periodic boundary conditions in all directions, 
on nanocrystalline Cu at 300 K. The normal stresses in the 
x and y directions perpendicular to the tensile direction z are 
maintained at zero, while the dimensions along these two di- 
rections are subjected to free adjustment. The system exists in 
the isothermal-isostress NcrT ensemble 12211 . Atoms interact 
with one another via the embedded atom model (EAM) Mishin 
potential 123l 12411 . The strain-rate in the tensile direction is 



3. Validation of the Nanocrystalline Model Composed of 
Randomly Packed Uniform Grains 

3.1. Geometrical Properties of the Grain Core and Grain 
Boundary Components 

3.1.1. Identifying Atoms within the Grain Boundaries 

In general, atoms located between the interface of two grains 
are treated as grain boundary atoms. To quantitatively mea- 
sure the volume fractions of the grain boundary and the grain 
core components, we use two methods to sort out grain bound- 
ary atoms. The first method is calculating coordination num- 
ber (CN). Here, atom i is a neighbour of atom j if the dis- 
tance r,j between them is smaller than r curo ff, which is de- 
fined by the value of the first minimum of the pair distribu- 
tion function of the system. For Cu at room temperature in our 
study, r cuto ff = 0.309 nm. The total number of neighbours of 
an atom is called the atom's coordination number. Atoms of 
CN < 12 are selected as grain boundary atoms and the rest are 
core atoms. The second method is Common Neighbour Anal- 
ysis (CNA). In this method, we apply the same r cuto ff for the 
conventional CNA calculation using LAMMPS SHEES] and 



also try the adaptive CNA (a-CNA) calculation using OVITO 
i28l 12911 . Fee and hexagonal close packed (hep) atoms are cate- 
gorized as grain core atoms. Both CNA and a-CNA calculations 
yield almost the same amounts of grain boundary and grain core 
atoms. 

3.1.2. Measuring the Average Thickness of the Grain Bound- 
aries 

Our Voronoi-like protocol generates grains shaped like con- 
vex polyhedrons. We approximate these convex polyhedral 
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grains by spheres with an average diameter d and an average 
volume V satisfying the relation 



3 \2 n 



core r core 



(4) 



where d is understood as dcN {dcNA) if N core , the total number 
of the grain core atoms, is obtained by the CN (CNA) calcu- 
lation, n is the total number of Voronoi grains, and V core is the 
average volume occupied by a single atom in an fee or hep unit 
cell, which is 0.0118 nm 3 for Cu at 300 K. Throughout our 
investigation, dcN - dcNA is always smaller than 0.4 nm. 
We can define the volume fraction of the grain core as 



nV/V„ 



(5) 



where V system is the volume of the simulation supercell. More- 
over, we assume that each approximated spherical grain is cov- 
ered with a grain boundary shell of average thickness dcB- The 
overall grain boundary thickness of the system is thus 2dcs, be- 
cause grain boundary is the space between two adjacent grains 
(see figure HJb)). Following the above assumptions, <p core can 
also be given by the expression 



(6) 



d + 2d C B 



Varying the size of the simulation supercell, we can create a set 
of (d, (p cor e), where d and <p core can be calculated using equa- 
tions and (0 separately. Both CN and CNA calculations 
give similar (p core when d > 10 nm. As a result, when greater 
than this system size, dcN and dcNA basically can be used inter- 
changeably and labeled simply as d. By fitting the set to equa- 
tion (O, we can estimate the grain boundary thickness in our 
system as shown in figure|5] We found that the total fitted grain 
boundary thickness of the system (2dcB) is about 0.39 nm (CN 
calculation) and 0.64 nm (CNA calculation), which are close to 
the values 0.5 nm) reported in other experimental and sim- 
ulation studies H3, [31, 32 , 33 ] . We also directly calculated 
the grain boundary thickness by dividing the total volume of 
the grain boundary component by the total interfacial area be- 
tween grains calculated by Voro++ i20ll . The thickness is about 
0.47 nm (CN calculation) and 0.53 nm (CNA calculation), as the 
grain size d is 25 nm, where E approaches the coarse-grained 
limit. 

3.1.3. Geometrical Inhomogeneity of the Grain Boundary 
Component 

In figure [6] we show the grain boundary structure. The 
grain boundary is about two to three atoms thick, translating 
to » 0.5 nm for Cu at 300 K. We can also observe some ge- 
ometrical inhomogeneity in the grain boundary region. This is 
an intrinsic feature of the grain boundary, enhanced by thermal 
fluctuations when the temperature is above zero. To analyze the 
sparsity of grain boundary atoms due to this inhomogeneity, we 
measured the average volume surrounding one atom in the grain 
boundary, V CB . The value is 0.01442 + 0.001 nm 3 (CN calcula- 
tion) and 0.01234 + 0.0002 nm 3 (CNA calculation), which are 




d (nm) 

Figure 5: Volume fraction of the grain core <f> cor€ as a function of d (d = 
dcN :black, d = rfcA^blue) of nanocrystalline configurations generated by the 
new Voronoi-like algorithm (solid circles), den - dcNA < 0.4 nm through the 
whole range of d. We use n = 8 for d > 12 nm, n = 27 for 4 nm < d < 12 ran, 
and n = 512 for d < 4 nm to keep the normalized standard deviation 
c s (dcN)/dcN ^ 0.1. Solid lines are least square fits of grain boundary thick- 
ness doB to Eq. |6}. The estimated overall grain boundary thickness {2dcB) is 
about 0.39 nm (CN calculation) and 0.64 nm (CNA calculation). Because dcB 
is near a constant, ty core — » 1 at the coarse-grained limit (d — > 00). 



22% and 5% larger than V core , respectively. Moreover, we cal- 
culated the average CN of the whole system (CNsy Stem ) and of 
the grain boundary (CNqb)- For dcN = 12.2 nm, we find they 
are 11.9 and 11.1, respectively, which agrees with a recent ex- 
perimental study with a similar d (d = 13 nm, CN S y Stem = 1 1.9, 
CNcb = 11.4) Q. 

3.2. Overall Elastic Response 

The overall E(d) is shown in figure [7] against experimental 
data EE HSU, and another MD simulation 11311 data. 
The experimental data are obtained within a temperature range 
of 300 K to 450 K, however. To compare our simulation data 
with the experimental ones at different temperatures meaning- 
fully, we measured the temperature-dependence of E in our MD 
tensile tests. The obtained value was -62 ± 7 MPa/K, which 
is almost the same as what had been reported in a previous MD 
study (-60+18 MPa/K) I3, but larger than one experimental 
value (-40 MPa[K) found in the same material with a grain 
size of 200 nm 13811 . For a temperature change of 150 K in sim- 
ulation, Young's modulus will vary by 10 GPa, which is also 
the size of the error bar that should be applied to the simulation 
data when we make a raw comparison. Similarly, we can calcu- 
late the projected values of the experimental data at 300 K if the 
temperature-dependence of -40 MPa/K is taken into account. 
The fitted E(d) agrees with experimental data extremely well, 
and approaches the coarse-grained limit 130 GPa lf39h when d 
is larger than 25 nm. 



3.3. Flow Stress of Plastic Deformation 

Using the uniaxial tensile test at 300 K with a strain-rate 
5 x 10 8 s _1 , we also measured the flow stress o> taken as the av- 
erage stress (cr) in the strain (e) interval between 7% and 10% 
as a function of d. In figure [8] we compare our simulation data 
of uniform grains with another simulation study, where a set of 
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Figure 6: 3D snapshot of the grain boundary structure of a system of d a 2 1 nm, 
where (dcN,<f> c ore) — (20.7 ran, 0.93) and (dcNA,<f>core) = (20.6 ran, 0,92). The 
dimension in the direction perpendicular to the figure plane is about one atom 
thick. We make all atoms of CN =12 invisible using Atomeye [3]. The inset 
is a blowup of the selected region. The grain boundary measures about two to 
three atoms thick with geometrical inhomogeneity. The configuration shown 
here may resemble an ordered packing visually, because the system contains 
only eight grains (n = 8) in 3D, only four of which are shown on the 2D local 
profile. Refer to figure^Ja) for its complete 3D Voronoi tessellation. 
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Figure 8: Flow stress o~p as a function of uniform grain size d compared with 
another simulation study using nonuniform Voronoi grains (black) [HJ]. We use 
larger n (8: blue, 27: green, 512: red) for smaller d to prevent premature grain 
growth during relaxations and tensile tests due to small grain size. Both studies 
give similar plastic strength and Hall-Petch-like behavior. 

Voronoi seeds are randomly chosen to create nonuniform grains 
1 1 Oil - Both studies show a transition from regular Hall-Petch 
to inverse Hall-Petch behavior around grain size d ~ 15 nm. 
Again, using uniform grains does not introduce extra weakness 
to our system due to monodispersity of the grain size distribu- 
tion. 

We also verified the isotropy of the system. For n = 8 
(smallest number of grains used in this study) grains of d « 
25 nm, we measured the elastic moduli and flow stresses in 
three perpendicular directions; their variations were only 6% 
and 8%, respectively. For n =512 grains, the system con- 
tains n x (n - l)/2 = 130816 different grain boundaries, and 
we believe it behaves isotropically. To make sure there is min- 
imal finite size effect, we measured the elastic moduli and flow 
stresses of n = 8 and n =27 grains of d ~ 12 nm, and n = 27 
and n = 512 grains of d ~ 5 nm. The relative variations in 
the elastic moduli and flow stresses in either case are smaller 
than 12% and 3%, respectively. In summary, all of the above 
agreements in the geometric, elastic, and plastic properties, and 
isotropy and finite size tests indicate that our new Voronoi- 
like protocol generates reasonably nanocrystalline configura- 
tions resembling real structures closely. 



Figure 7: The overall Young's moduli E as a function of d. The black dots (orig- 
inal MD data), black (CN calculation, fitted), and blue (CNA calculation, fit- 
ted) lines are from this study; the cyan circle I13|l represents another simulation 
data; the yellow squares (331, red squares ]35l |36|| . orange square @J, and blue 
square 11611 show data from experimental studies; the dashed line represents 
the coarse-grained limit 13911 . Color bars associated with the experimental data 
show projected values of the moduli at 300 K by considering the temperature- 
dependence (-40 MPa/K) observed experimentally I3H1 . if the original data 
were not obtained at this temperature. Color bar associated with the cyan sim- 
ulation data shows the offset radius dcN , because grain boundary thickness is 
not taken into account in that study. 



4. Analysis of the Elastic Moduli of the Grain Core and 
Grain Boundary Components 

We can express the overall Young's modulus E in terms of 
(pore, E core and Eqb by rules of mixtures, where E core and Eqb 
are Young's moduli of the grain core component and the grain 
boundary component, respectively. To estimate E core and E<jb, 
we fitted E to the Voigt Reuss JUl, and Hill 01 models, 
representing the upper bound of E, the lower bound of E, and 
the arithmetic mean of the previous two models, respectively. 
The Hill model has been chosen for practical reasons 1I43I . For 
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Table 1: Young's moduli (E core , Egb, Egb I 'E core ) 



(Unit: GPa) 


Voigt 


Hill 


Reuss 


(2- 


norm) 










CN 


(138, -125, <0) 


(167, 10,0.06) 


(177,25,0.14) 




CNA 


(130, -8, <0) 


(158,20,0.13) 


(153,40,0.26) 


(oo — 


norm) 
CN 
CNA 


(136, -130, <0) 
(127, -7, <0) 


(124,26,0.21) 
(136,32,0.24) 


(145,32,0.22) 
(146, 40,0.27) 




Figure 9: Numerical fits with oo - norm (solid lines) of the overall Young's 
moduli E of <f> core (solid circles) determined by CN (black) and CNA (blue) 
calculations to the Reuss model. A hypothetical Reuss curve l4lh assuming 
Ecore = 130 GPa and Egb = 0.7E core (dashed line) is also shown to emphasize 
the disagreement between the hypothesis and the simulation data of this study. 



a given set of E and d> core , E core and Ecb become parameters to 
be determined by numerical fitting. We tried both the 2 - norm 
and oo - norm. The results are listed in Table Q] 

The Voigt model gives unphysical negative values of Eqb in 
all four cases. Under both definitions of grain boundary atoms, 
Hill and Reuss models give meaningful numbers worth further 
consideration. Theoretically, the Reuss model describes the ef- 
fective moduli of a solid suspension in a matrix of near zero 
shear modulus, a quantity proportional to the elastic modulus 
for homogeneous isotropic materials [14411 . In Table Q] we found 
that all ratios of EcBlE core are smaller than 0.27, a close sce- 
nario to an ideal Reuss example, although the ratio is nonzero. 
In other words, in our system, grain cores are embedded in a 
grain boundary matrix with a stiffness of less than about 1 /4 of 
that of the core part. This result reciprocally justifies the fact 
that among the three proposed models, the Reuss model best 
captures the elastic response of our system. Specifically, as- 
suming a linear constitutive relation, the Reuss model can be 
expressed as 



E = 



1 



/Ecore + (1 - 4>core)/EcB 



(7) 



The Reuss fits are shown in figure [9] We choose the results of 
oo - norm over 2 - norm, because the E core values are much 
closer to the coarse-grained limit, 130 GPa i39ll . 

Knowing cLgb, E core , and Eqb, we can rewrite £ as a function 
of only d by combining equations (O and (0. The fitted overall 




Figure 10: Probability distribution of the least square von Mises shear strain 
1431 of the grain core (dashed line) and the grain boundary (solid line) obtained 
by CN calculation under thermal fluctuation at 300 K. d « 11 nm, where 
(dcN,<pcore) = (11.1 nm, 0.89) and (dcNA^care) — (10.9 nm, 0.85) in this sys- 
tem. The results are obtained by comparing two snapshots separated by 10 ps. 
The grain boundary is subjected to more deformation even though there is no 
external loading. 



E(d) is shown in fig ure[7 ]compared with other experimental and 
simulation data OElsjlMlIlllEl. 

The greater temperature-dependence in our system and the 
elastic softness of the grain boundary component may be be- 
cause the volume fraction of the grain boundary increases 
with decreasing d, and the grain boundary is more thermally- 
sensitive than the grain core, a fact due to the geometrical inho- 
mogeneity in the grain boundary (see the inset of figure |6) lll3ll . 
We verify this by looking at atomic scale deformation due to 
thermal fluctuations of both components at 300 K. In figure [TOl 
we show that, locally, the grain boundary deforms more than 
the grain interior, while both components are subjected to an 
equal amount of thermal vibrations at the same temperature. 



5. Discussions and Conclusions 

In this study we investigated Young's moduli of the grain 
boundary and the grain core components of fully dense 
nanocrystalline copper with an average grain size d < 25 nm 
at 300 K. We used an isotropic nanocrystalline model that 
includes not only the most stable grain boundary, but also 
a general combination of varied grain boundaries generated 
by our new Voronoi-like algorithm. The seeds of the new 
Voronoi algorithm, used for creating monodisperse nanocrys- 
talline grains, are centers of randomly close packed identical 
spheres. The model composed of uniform Voronoi grains, al- 
though extremely ideal, faithfully reproduces the geometric, 
elastic, and plastic features of bulk nanocrystalline copper. We 
conducted MD uniaxial tensile tests to measure Young's modu- 
lus as a function of d. 

We found the following key results concerning the stiffness 
of nanocrystalline Cu at finite temperature: 1) The grain bound- 
ary component is less stiff than the grain interior; 2) At 300 K, 
Young's modulus of the grain boundary component is less than 
30% of that of the grain interior, not 70% as generally be- 
lieved; 3) The grain boundary component is more thermally- 
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sensitive compared with the grain core component, because the 
grain boundary shows more local deformation compared with 
the core region under the same thermal fluctuations. 

It is important to emphasize here the distinct simplicity of our 
approach to estimating Young's moduli at finite temperature. 
Starting with a nanocrystalline model containing near-spherical 
monodisperse grains with random crystal orientations, we ob- 
tain the grain boundary thickness and its elastic modulus by nu- 
merical fits, without making any assumptions about the elastic 
responses of the grain boundary. Due to unavoidable inaccu- 
racy of numerical fits, we focus on the ratio of Young's moduli 
in the grain boundary and the core components, instead of their 
absolute values. Through our analysis, we argue that the grain 
boundary is elastically much softer, with less than 30% of the 
stiffness of the grain interior, at the studied temperature. This 
striking finding could be explained by the fact that the some- 
what disordered grain boundary component is much more ther- 
mally sensitive, and it has been verified by observing the struc- 
ture change under thermal fluctuations (figure [10). Because the 
CN of the grain boundary component is only slightly smaller 
(< 7% for d « 12 nm, for example) than that of the whole sys- 
tem and the volume shared by a single grain boundary atom is 
merely 5% to 20% bigger than that of a core atom, it is impres- 
sive that a small disorder in structure can create a significant 
drop in stiffness. 

Studies assuming Young's modulus of the grain boundary to 
be 70% of the grain interior value usually attribute the observed 
reduction of the overall elastic modulus to porosity when the 
grain size is smaller than 20 nm JsJ. The influence of poros- 
ity on the elastic properties of nanocrystalline materials at finite 
temperature is likely overestimated and should be further exam- 
ined, since we have shown here that the grain boundary is not 
as stiff as previously believed. 

Our results also directly help developing more advanced con- 
tinuum models. For example, in recent quasi-continuum simu- 
lations, it has been shown that if the grain boundary has a very 
small elastic modulus (17% of the stiffness of the grain core), 
the shear stress field of the edge dislocation shrinks with grain 
size 114611 . Our study about Young's modulus in the grain bound- 
ary region sheds new light on this result. 
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